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Fig. 1. Buckling of the C11O2 plane. O atoms are displaced from 
the CUO2 plane. Arrows represent the spins at Cu sites. The 
hole at the O site interacts with the spin via Kondo coupling 
J K . 



We show that a Chern-Simons term for a gauge field describing a fluctuation of spins is induced 
by integrating out hole fields in the presence of spin-orbit coupling which originates from a 
buckling of the Cu02 plane. Through the Chern-Simons term, holes behave like skyrmion 
excitations in a spin system and become a superconducting state with d x 2_ y 2 symmetry after 
the antiferromagnetic long-range order is destroyed. 

KEYWORDS: high-T c superconductivity, mechanism of d-wave superconductivity, Chern-Simons term, skyrmion 
excitation, buckling of CuC>2 plane, spin-orbit coupling 



Since the discovery of high-T c superconductivity in 
cuprateSjEP much experimental and theoretical effort has 
been invested to clarify its mechanism of superconductiv- 
ity. Results of experimental studies indicate that the fol- 
lowing properties are essential features. First, the CuC>2 
layered structure is intrinsic to superconductivity and 
both the undoped and carrier-doped CuC>2 planes are 
characterized as, two-dimensional systems on the basis of 
their magnetic,!^ transport^ and optical properties.cP 
Second, superconductivity occurs in a disordered spin 
background.EF For the undoped case, the system is a 
Mott insulator and spins at Cu sites show antiferromag- 
netic long-range order below the Neel temperature T/\r. 
Upon doping, TV decreases to zero and spin-glass be- 
havior is observed.© Superconductivity emerges upon 
further doping. Apparently, disorder in the spin system 
is introduced by doped holes. Third, the Cooper pair is 
spin-singlet and has d x ?_ y i symmetry.© 

In addition to these properties, the occurrence of su- 
perconductivity appears to be closely related to the 
structure of the Cu02 plane. Near the overdoped re- 
gion of La2_ x Sr x Cu02 the disappearance of supercon- 
ductivity was observed at an ortharhombic to tetragonal 
structural phase transition point u Moreover, the tran- 
sition temperature is closely related to the buckling 
of the CuC>2 plane.tJ Among the effects on the conduc- 
tion electron system accompanied by a biiddipg of the 
Cu02 plane, there is spin-orbit coupling £3' BO Spin- 
orbit coupling can have an important effect on conduc- 
tion electrons through the Bern/ phase induced by the 
background spin configuration.EiJ 

In this Letter, we propose a mechanism of d-wave su- 
perconductivity in a disordered spin background based 
on a two-dimensional model with spin-orbit coupling. 
We assume for the spin-orbit coupling term that it is 
induced by the buckling of the Cu02 plane, as shown in 
Fig. 1. We show that a Chern-Simons term for a gauge 
field, which describes the fluctuation of spins at Cu sites, 
is induced by integrating out the hole fields. Through 
this Chern-Simons term, holes behave like skyrmion ex- 



citationst3'li2Hi2Hiil't2Hi^ for the spin system. When the 
antiferromagnetic long-range order is destroyed by these 
skyrmion excitations, the Chern-Simons term leads to 
Cooper pairing of holes. We show that the pairing state 
is spin-singlet with d x 2_ y 2 symmetry usirig.a transforma- 
tion to the previously considered modcl.EI 

Our model is described by the following Hamiltonian: 



H 



-*<>£ 

(i,3) 



H.c 



-j k £.. 



where c. 



( <4 T cjj J and Ci = ( c lT Cii 



(1) 
) T de- 



note a creation and an annihilation operator for holes in 
a spinor representation. In the first term, the summa- 
tion is taken over the nearest-neighbor sites. The sec- 
ond term denotes Kondo coupling between the spin.-cjf 
holes: = \c\ctc^ and the spin at Cu sites: SjE3 



3 r 2 t j 

The Hamiltonian 



is based on a model which distin- 
guishes electrons at the Cu site and holes at the O site. 
We implicitly exclude the double occupancy of holes at 
O sites because there is a strong on-site Coulomb re- 
pulsion. A similar model without i7 so is proposed in 
refs. ^l], ^2] and ^3|. Note that if we take the limit 
Jk — > 00, eq. (Oli jssdthout H so corresponds to the so- 



called t-J model 



For spin orbit coupling, we as- 
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sume the following formal!}!! 

H *o = i J2H c l x(a) ■ a c i+™<* + H.c, (2) 



j a=x,y 



where A^ — (x^, ^y*j , e Q is a unit vector along the 
ct-axis, and a is the lattice constant. Spin-orbit coupling 
(||) is produced by a buckling of the CuC>2 plane. An 
example is shown in Fig. 1. In the presence of spin- 
orbit coupling (||), there is generally a Dzyaloshinskii- 
Moriya-type interaction for the spin system. However, 
we ignore this term because it does not play an impor- 
tant role in terms of the mechanism of interest. In order 
to describe the spin system we introduce the Schwinger 
bosons© Sj = \z\crz^ with z\ = { z\+ z\ t ^ and 



l t 

zj = ( zji Zji ) . Here Zj a are boson fields and obey 
the constraint z^Ztf + z ]i z ji — 1> since S — 1/2. 

We use the path-integral formulation to describe the 
system. The action is given by 

' (3) 

The effect of the spin fluctuation on holes can be intro- 
duced by performing a series of SU(2) transformations. 
Let us focus on an i — j bond along the a-axis where 
i and j are nearest-neighbor sites. We assume that the 
j-site belongs to the A-sublattice and the i-site belongs 
to the B-sublattice. The hopping term is given by 



^(A (a) ) = ^^(A (a) )c 
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(4) 



(ah 



iX 



(a) 



<7 ) , with oo 



where y(A , 

the unit matrix in spin space and X a = |A^*'|. Since 
X a =/= 0, there is a twist in spin space. Then, first we per- 
form a transformation to diagonalize the hopping term: 



X 



(A^)cj, and Cj 



c,;. By this transformation, 
up toO((A Q /t ) 2 ) 



X 2 a ciCj 



we obtain hij(X {a >) = - V^o 
in the phase factor. Second, we perform a transfor- 
mation which diagonalizes the Kondo coupling term 

cj — > UjCj, Ci — > Ci (icfyUi), where Ui = ( Z ' T 

/ ~zt\ ~z~i i \ 

and Ut = ■ By this transformation, 

V ~ z n z n J 

SU(2) gauge fields wMch^describe the fluctuation of the 
spins are introducedlH§EZl ) UJJj = \UiUj\ exp (— iaAji) . 
The use of this transformation is based on the assump- 
tion that Jk has a higher energy scale than any other 
parameter in the Hamiltonian and the antiferromagnetic 
correlation between spins is retained and its correlation 
length is much larger than the lattice constant a. In this 
transformation, the sign of Kondo coupling and the time 
component of the gauge field: At = —iUdtU , for holes at 
the B-sublattice is changed since there is an extra factor, 
iuy for the B-sublattice. For the amplitude fluctuation of 
the spins; \UiUj\, we take its mean value rj = (\UiUj\), 
because it is a relatively higher energy mode than the 
phase fluctuation of the spins. Next, we perform the in- 
verse transformation of the first transformation and we 



obtain 



x exp 



1 \(<*) A 

A v ' ■ a — iclA 



U 



3i 



3 ' 



(5) 



up to 0((A Q /io) 2 ) in the exponent. 

Compared with the model considered in ref. the 
derivation of the Chcrn-Simons term is complicated be- 
cause we need to distinguish between the A- and B- 
sublattices. However, there is a transformation which 
connects the model to that in ref. BO, that is, 



3 ' 



Ci — > Ci(-ia y ). 



(6) 



After this transformation is performed, an additional 
sign change occurs for both the Kondo coupling term 
and the time component of the gauge field. However, as 
we will discuss later, transformation (^) changes the sym- 
metry of the pairing matrix. Therefore, we must perform 
the inverse transformation to obtain the physical state of 
the original system. In the following, we call the system 
obtained after transformation (^) the "F-system" and 
the original system the "AF-system." In the F-system, 
the action for the holes is given by 

S h = f dtJ2c,(t)G^({k ~ ^})c,-(*), (7) 
J j 

where fc p is defined by e laka Cj — Cj+ a & a , for fi = x,y, 
and kt = —idf The inverse of the Green's function is 
given by 

G- 1 {{k^}) = {k + 2t oV cosfc Q )a -g(k) • <x, (8) 

a—x,y 

where 



g(k)= 2r/ Y, Ai Q) sinfc Q ,277 £ A^sin^, 



Jk 



a—x,y 



(9) 

Note that we cannot take the limit Jk — * oo in the 
presence of spin-orbit coupling (|J). If we take the limit 
Jk — > oo, then the spin of holes is projected in the di- 
rection antiparallel to the spin at Cu sites. However, 
the hopping process always involves the opposite spin of 
holes as long as A Q ^ 0. 

In order to calculate the Chern-Simons term, we take a 
continuum limit. (The condition of taking this limit will 
be discussed, later.) The induced Chern-Simons term is 
given byE§E3'l23) 



S cs = -^JdtJ d 2 vA x t {d x A x y - d y A%) 



where 



— /^'(^gfklx^gtk)) 

47T 



|g(k) 



(10) 



(11) 



Note that only the Abelian Chern-Simons term appears 
in eq. ( |l0| ) because the SU(2) gauge field A^ is reduced 
to the AbeliaruChern-Simons term upon using the curl- 
free condition^ We retain Af, because it describes the 
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staggered spin fluctuation. From eq. (|l l|) , we obtain 
1 



x sgn (J K A) , 



(12) 



!3D 



where A = Ax Ay — Ay Ax , and we have used the con- 
tinuum form of (j|) : g(k) = {2r)J2 a x x a) ko^2n^2 a A y a) k a , 
— Jjf/4). In contrast to the anyon system,tifEf the va 
of 8 = ±1/2 does not alter the statistics of particles. 
We can extend the above calculation to a finite tem- 
perature.EJ However, if we concentrate on the region 
ksT <C Jk, we can neglect finite temperature correc- 
tions. Since spin-orbit coupling term (j^) involves a pro- 
cess of hopping between different d-orbitals at the same 
site, the external electromagnetic gauge field does 
not couple to it. Therefore, there is no Chern-Simons 
term for A° xt . In deriving eq. (|Io|), we have taken the 
continuum limit for the gauge field A^. Since the length 
scale of the gauge field is given by u/A sw , where A sw and 
v denote the gap and the velocity of the spin wave mode 
respectively, the condition of taking the continuum limit 
is A sw /(v/a) < A a /to, which can be seen from eq. (5). 
For A a , a rough estimation gives A a ~ 2 meViliJ We 
assume that this condition is satisfied in the underdoped 
region because there A sw may be very small and v/a-is 
close to the value of the undoped case ~ 200 meV.lia) 
Although the presence of spin-orbit couping is essential 
for the derivation of the Chern-Simons term, it has no 
importance for other physical processes. Therefore, we 
can neglect it in the following discussion. 

Now we discuss the effect of the Chern-Simons term. 
We take 



S — Sh + S, 



cs 



S, 



spin j 



(13) 



for the effective action. The last term is the action for the 
spin system and is given by the CP\ modekES 1 5 sp i n = 



(2/ff)/d 3 a:£„ \{d, - iAfi z„\ z + (A s 2 w /« 2 ) \z a f 



By 



integrating out the Af, we obtain the relationship be- 
tween the spin density of the hole and the "magnetic" 
field B(r,i) = d x A x y {v,t) - d y A%{r,t). If we take the 
x-axis as the quantization axis for the spin, we obtain 
Y, a SaPa{r,t) = j^B(r,t), where s a = 1 for a =\ and 
s a = —1 for a =[. In the AF-system, this relation in- 
volves the isospin index, that is, 



E 



s a s T p a (r,t) 



2tt 



B(r,t). 



(14) 



Here s T = 1(— 1) for r belongs to the A(B)-sublattice. 
Therefore, T(j)-spin at the A-sublattice induces a (anti- 
)skyrmion excitation in the localized spin system and 
!(|)-spin at the B-sublattice induces a (anti-)skyrmion 
excitation in the localized spin system. Since the 
skyrmion and anti-skyrmion excitations introduce dis- 
order into the localized spin systemEZP and the number 
of them is the same as that of holes, disorder in the 
spin system increases upon doping. If the hole density is 
sufficiently small that the magnetic long-range order is 
preserved, then the Meissner effect occurs for the gauge 
field A^ and holes are pinned because skyrmions break 
the translational invariance, that is, the system is an in- 
sulator. 



After the magnetic lonfeKange order is destroyed by 
the skyrmion excitations,t3 , tH' the Chern-Simons term 
becomes dominant in the long wavelength and low- 
energy physics. For the holes, it leads to a pairing state. 
Coupling between the hole current and the gauge field 
A x , in the F-system is given by 



S 3-A = 



E 



J dt J d 2 r s a j a (r,t)-A x (r,t), 



(15) 



where j ff (r,t) is the hole current for c-spin. Since eq. 
( |l5| ) describes minimal coupling between the hole current 
and the gauge field A. x , it gives rise to a Lorentz force. 
Such a Lorentz force is induced between holes passing 
each other. Therefore, it leads to a chiral pairing state. 
The chirality is determined by the sign of 9. From its 
pairing mechanism, the possibility of the s-wave pairing 
state is excluded. 

Now we investigate the pairing state of the AF-system 
through the F-system. Before doing that, we must know 
the relationship of the pairing matrix between them. We 
assume that i and j are nearest neighbor sites. If we 
take Af ■ = (c^Cji — CjjCjf) for the spin-singlet pairing 
order parameter for the AF-system, then after perform- 
ing transformation (||) we obtain A|. — > (c^Cj^ +Cjj.Cjj.). 
Therefore, the spin-singlet pairing state is transformed 
into the spin-triplet pairing state and vice versa. In 
k-space, holes at the A-sublattice are described by the 
fields aj^ = f c k CT + c k+Qo-) anc ^ n °l es at the B- 

sublattice are described by = f Cj C(T — C k+Q CT ) > 
where Q = (jr/a,Tr/a). Here, we assume that the A- 
sublattice is the set of (even, even) and (odd, odd) and 
the B-sublattice is the set of (even, odd) and (odd, even) . 
The pairing matrix may be given by 



A 



E^- 



kv 2 a k' ffl MF, 



(16) 



where Yfcifk + /k+o) = Ek /k- Here > we do not need 
the explicit form of , because we solve the gap equa- 
tion through that of the F-system. Transformation (||) in 
k-space is given by (3^ — > iPyP u . By performing this 



transformation, the gap equation (16) is transformed into 



^ 1CT2 =E V kk((^-k 



For the singlet pairing case A^ 
reduced to 



n 



-A 



IT' 



IF- 



eq. 



A 



n 



A 



(i)k , A ( 2 )k 



TT 



TT 



/2, 



where 



A TT k = E V kk( c -kT c kT) F ' 



(17) 



(18) 



(19) 



A (f = -EV k F k ,(c_ k , +QT c k , +QT ) F . (20) 
k' 

The minus sign in eq. (20) originates from the 
sign change in the kinetic term in G _1 (fco, k + Q). 
Although the vector g(k) changes as g (k + Q) = 



d 
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diag(— 1, — 1, l)g (k), the sign of 9 does not change. In 
the continuum approximation, eqs. (19) and (20) are 
reduced to 

4-7ri k x k' A k 

A k = ±E-TikTW2^' (21) 
k' 



where = 



|A k l 2 , 



at T = 



analysis of ref. 
tion is given by A*i?^ 



0. Following the 
we can solve eq. pi]). The solu- 
A fe exp (±24%), and A^ )k = 
Afeexp hf2i£6] s ), where A/, is a function of k = |k| and 
^k = arc t an ^i//fca:- Here, 21 is the relative angular mo- 
mentum of the Cooper pair and is not equal to zero. The 
gap of superconductivity A^ F is of the order of Fermi 
energy tp and the cohrence length of superconductivity 
&3C is given by £ S c/a ~ 4%/ | Afc ■ Here > x is the hole 
concentration. Since Afc p is of the order of ep, this value 
may be smaller than the average distance between holes: 
~ a/*Jx. From eq. (|i"8|), we obtain 



A T1-- A 1T= A ^ 0S ( 2 %) 



(22) 



Since the smallest I is realized in the ground state, we 
set I = 1. In this case, eq. ( p2| ) describes the d x 2_ y i 
pairing state because cos (20^ J = (k x — k y )/k 2 . For 
the triplet pairing case, we find that A^ = A^ = 



( 



,(i)k 



Al 2)kN l /2, and A k 



= A k T = 0. 



However, such 



a pairing state is not stable in the bulk of the system be- 
cause the d- vectored satisfies d^. || e y , that is, the spins 
of Cooper pairs lie in the plane perpendicular to the y- 
axis. As a result, the pairing state has spin-singlet and 
d x 2_ y 2 symmetry. 

There is also some contribution to the pairing mech- 
anism from other spin fluctuations, which may be char- 
acterized by the Maxwell term: ~ — j (d^A* — d v A^) 
in the gauge field description because coupling to the 
spin system is mediated by the gauge field A*. Mean- 
while, our spin fluctuation is characterized by the Chern- 
Simons term. Since there is an extra derivative for the 
former compared with the latter, our mechanism may 
be more dominant in the long wavelength and the low- 
energy limit than other spin fluctuation mechanisms. 
Moreover, the Chern-Simons term only exists in the 2 + 1 
dimension. Therefore, our spin fluctuation is unique to 
the 2 + 1 dimension. In contrast, the Maxwell term exists 
in any dimension. 

For the apph^aliaa to the orthorhombic phase of 
La2_ 2; Sr 2 ;Cu04,E3Ot3) we require one more transfor- 
mation after transformation (||), that is, C( 0( jd.cvcn) - » 
icr z c (odd,cvcn), and C( ovcn odd ) — > —i<7zC(cvcn.odd)- Also in 
this case, the value of 9 is given by eq. JI4 ) and the 
pairing state has spin-singlet and d x 2_ y 2 symmetry. 

In summary, we have studied a model of the Cu02 
plane with buckling and have shown that the Chern- 
Simons term for the gauge field, which describes the 
fluctuation of the spin system, is induced. Through 
this Chern-Simons term, the doped hole behaves like a 
skyrmion or anti-skyrmion excitation depending on its 
spin or isospin, that is, whether it resides on the A- 
sublattice or B-sublattice. After the antiferromagnetic 



long-range order is destroyed by the skyrmion excita- 
tions, the Chern-Simons term becomes dominant for long 
wavelength and low-energy physics and leads to the spin- 
singlet d x 2_ y 2 superconducting state. 
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